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STRUCTURE  AND  DYNAMICS  OF  CrCHjO)*  CLUSTERS.  THE  EFFECT  OF  THE 
POLARIZABILITY  AND  THE  CHARGE  OF  THE  ION. 

Lalith  Perera  and  Max  L  Berkowitz 

Department  of  Chemistry,  University  of  North  Carolina,  .  Chapel  Hill, 

North  Carolina,  27599 


ABSTRACT: 

The  effect  of  the  polarizability  and  the  sign  of  the  ionic  charge  were  studied  in 
CTCHjOa,  clusters  using  molecular  dynamics  computer  simulation  technique.  From  our 
simulations  we  concluded  that  the  reduction  in  the  ionic  polarizability  did  not  significantly 
change  the  structure  and  dynamics  of  Cl'CHjO^  cluster,  but  the  inversion  of  the  sign  of  the 
ionic  charge  produced  a  large  effect  The  energetic  considerations  helped  us  to  understand 
why  Cl*  is  located  on  the  surface  of  the  cluster.  By  being  on  the  surface  the  anion  permits 
the  creation  of  the  hydrogen  bonded  network  between  water  molecules  and  that  lowers  the 
total  energy  of  the  cluster. 

Simulations  with  the  inverted  sign  of  the  ionic  charge  correspond  to  that  with  a 
hypothetical  "Cl*"  ion  which  is  sirai!**  in  size  and  polarizability  to  a  Ca*  ion.  The  dynamical 
structures  and  the  quenched  structures  of  *C1+"(HjO)jo  clusters  are  compared  with  the 
idealized  structure  of  Ca* (I^O)®  cluster  proposed  recently  (A  Se linger  and  A  W.  Castleman, 
Jr.,  J.  Phys.  Chem.  95,  6442  (1991)). 
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I  INTRODUCTION 

Recent  experimental  and  theoretical  studies  performed  on  aqueous  clusters  display  a 
rich  dynamical  and  structural  behavior  that  characterises  these  systems  [1-111.  Some 
experiments  indicate  that  water  structures  of  well-defined  geometry  may  enclose  ions,  such 
as  NH4*  [4],  HsO*  [5,6],  Off  [7]  and  Cs*  [8].  Of  special  interest  are  clusters  where  20  water 
molecules  are  encaging  the  ion.  It  was  proposed  that  the  enhanced  stability  of  these  clusters 
is  due  to  the  formation  of  a  pentagonal  dodecahedron  which  surrounds  the  ion  [6,8].  In 
contrast,  other  experiments  indicate  that  water  structure  in  a  cluster  solvating  anions  such 
as  Cl*  may  not  be  that  regular  [9].  Moreover,  some  of  the  cluster  structures  in  this  case  are 
rather  unexpected.  For  example,  in  our  recent  molecular  dynamics  study  of  G*(H,0)14  cluster 
we  obtained  a  rather  interesting  result:  in  this  cluster  the  anion  is  partially  solvated,  it  is 
located  on  the  surface  of  the  duster  [10],  Similar  condusions  on  the  character  of  solvation 
were  also  reached  in  the  studies  performed  on  electron  solvated  in  water  dusters  [11,12]. 

The  immediate  question  that  comes  to  one’s  mind  is  why  the  Cl*  ion  is  situated  on  a 
surface  of  the  duster?  Is  that  due  to  the  large  polarizability  of  the  Cl*  ion,  or  is  it  due  to  the 
sign  of  the  ionic  charge?  To  answer  these  questions  we  performed  three  molecular  dynamics 
simulations.  In  all  three  simulations  the  system  was  composed  of  20  water  molecules  and  one 
ion.  In  the  first  simulation  we  considered  the  dynamics  of  20  water  molecules  and  one  Cl*  ion; 
in  the  second  simulation  we  studied  the  structure  and  dynamics  of  a  duster  of  20  water 
molecules  and  one  "Cl*"  ion,  i.  e.  in  this  system  the  ion/water  interaction  was  described  by 
the  same  set  of  parameters  as  in  the  previous  simulation,  except  the  charge  on  the  Cl*  ion  was 
changed  from  -1  to  + 1.  By  doing  this,  we  were  able  to  investigate  how  the  sign  of  the  ionic 
charge  affects  the  solvation  of  the  ion.  It  is  appropriate  to  mention  here  that  the  "Cl*"  ion 
is  very  similar  to  the  Cs*  ion,  since  the  sizes  and  polarizabilities  of  the  two  ions  are  very 
similar.  (The  ionic  radius  of  Cl*  is  L81  A,  while  the  ionic  radius  of  Cs*  is  1.7  A.  The 
polarizabilities  of  Cl*  and  Cs*  are  3.25  A3  and  2.56  A3  respectively). 

In  order  to  study  how  the  large  polarizability  of  the  Cl*  ion  influences  the  dynamics 
and  the  structure  of  the  duster,  we  reduced  the  value  of  Cl*  polarizability  in  the  third 
simulation.  In  what  follows,  the  Cl*  ion  with  a  reduced  polarizability  is  referred  as  Cl*^  ion. 

H  METHODS 

In  our  previous  work  we  have  shown  that  the  indusion  of  many  body  interaction 
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effects  into  the  molecular  dynamics  simulation  of  ion/water  dusters  is  essential  in  order  to 
reproduce  the  correct  energetics  of  the  dusters  [10].  Therefore  in  these  simulations  we  again 
used  the  SPCE/POL  potential  model  which  indudes  many-body  effects  [13,14].  In  this  model 
the  water  molecule  is  assumed  to  have  three  charges  placed  at  the  positions  of  two  hydrogens 
and  the  oxygen.  The  Lennard-Jones  center  is  placed  on  the  oxygen  site.  The  OH  distance 
in  the  model  is  1 A  and  HOH  angle  is  109.47*.  Three  polarizable  centers  are  placed  at  each 
atomic  site  in  water.  The  ion  in  this  model  is  considered  to  be  a  point  charge  placed  at  a 
Lennard-Jones  center.  Additionally  a  polarization  center  is  also  placed  at  the  same  location. 

The  total  energy  of  the  system  is  given  by  the  following  expression: 

um  -  UsU'j.U^U^  (1) 


where 

-  YtMj/ra  (2) 


is  the  sum  of  all  Coulomb  interactions  between  all  i  and  j  charge  sites  and 


uu  -  E 

€<• 


(3) 


is  the  sum  of  the  Lennard-Jones  interactions  between  a  and  S  Lennard-Jones  sites.  The 
polarization  energy  is  obtained  from  the  expression  [15] 


(4) 


where  a,  is  the  polarizability,  /i,  is  the  induced  dipole  at  the  polarizable  i-th  center  and  r0  is 
the  distance  vector  from  site  i  to  site  j.  The  values  of  the  induced  dipoles  at  the  polarizable 
centers  are  calculated  self-consistently  using  relationships  (5)  and  (6) 

(5) 
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Ei  -  Ef-Y.TM  (6) 

where  E,  is  the  total  electrostatic  field  at  the  i-th  center  and  E,°  is  the  field  due  to  permanent 
charges.  The  dipole  tensor  T8  in  eq  (6)  is 

-  l/rJU-Zr^/r/]  (7) 

The  potential  energy  due  to  the  three  body  exchange  repulsion  of  the  ion-water  trimer  is 
expressed  in  terms  of  ion-oxygen  distances  ra  and  r^  and  the  oxygen-oxygen  distance  r^; 

^3-Mr  “  Sexp< ~ )exp( - flri4)exp( -yrJh)  (8) 

where  B,  0  and  y  are  parameters.  The  potential  parameters  of  the  SPCE/POL  model  are 
given  in  Table  I. 

Molecular  Dynamics. 

We  used  molecular  dynamics  computer  simulations  to  investigate  the  structure  and 
dynamics  of  small  clusters  made  up  of  one  ion  and  20  water  molecules.  Equations  of  motion 
were  solved  using  the  Verlet  algorithm  [16].  The  geometries  of  the  water  molecules  at  their 
equilibrium  values  were  constrained  through  the  Shake  procedure  [17].  To  solve  self- 
consistent  equations  (6)  and  (6)  an  iterative  approach  was  used.  The  iterations  were 
continued  until  the  root  mean  square  of  the  difference  in  the  induced  dipoles  between 
successive  iterations  was  less  than  0.001  D/atom.  The  self-consistency  was  usualy  achieved 
in  seven  steps  of  iteration. 

Each  cluster  was  initially  equilibrated  for  more  than  20  ps  followed  by  500  ps  of  the 
production  run.  The  time  step  of  the  trajectory  has  been  kept  at  1  fs  and  the  total  linear  and 
angular  momenta  of  the  system  were  removed.  During  the  trajectory  calculations  the  energy 
was  conserved  within  5-th  significant  digit 

All  the  trajectory  calculations  were  done  in  the  microcanonical  ensemble  at 
temperatures  around  200  K  The  observed  temperature  fluctuations  were  around  30  K.  The 
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temperature  of  the  cluster  was  chosen  to  be  high  enough  to  eliminate  possible  trapping  of  the 
cluster  in  a  certain  configuration  for  a  very  long  time.  We  attempted  to  perform  some 
molecular  dynamics  simulations  of  our  clusters  at  300  K,  but  observed  the  evaporation  of 
water  molecules  at  this  temperature. 

m  RESULTS  and  DISCUSSIONS 

Since  the  net  translational  and  rotational  motions  of  the  center  of  mam  were  absent 
in  our  simulations,  we  used  the  distance  from  the  ion  to  the  origin  aa  an  index  of  the  cluster 
solvation.  If  a  duster  is  completely  solvated,  this  distance,  as  well  as  the  distance  to  the 
center  of  mass  of  the  water  molecules  should  approach  zero.  The  two  distances,  shown  in 
Figures  la,  lb  and  lc  also  provide  information  on  whether  the  duster  is  symmetric  when  the 
solvation  takes  place.  From  comparison  of  Figures  la,  lb  and  lc  we  can  see  that  the  most 
solvated  ion  is  ”C1*". 

To  gain  more  insight  into  the  structure  of  the  dusters  we  calculated  the  ion-water 
radial  distribution  functions  (rdf),  which  are  represented  in  Figures  2  and  3.  Figure  2  displays 
the  ion-oxygen  rdf  and  Figure  3  displays  the  ion-hydrogen  rdf.  For  Cl'/water  duster  the  first 
maximum  in  ion-oxygen  rdf  is  at  r«3.06  A;  the  coordination  number  of  the  ion  is  4.5. 
However,  for  Cl'  ion  solvated  in  bulk  water  at  300  K  these  values  are  correspondingly  3.2  A 
and  6.1  [14],  The  reduction  of  the  polarizability  of  Cl'  ion  produces  a  shift  in  the  position  of 
the  first  maximum  of  ion-oxygen  rdf  towards  a  larger  value  of  the  distance  and  an  increase 
of  the  coordination  number  to  the  value  of  5.5.  hi  the  case  of  Cl*  ion  the  radial  distribution 
functions  clearly  show  that  the  solvation  shell  consists  of  a  primary  shell  and  a  broad 
secondary  shell.  This  is  not  so  in  the  case  of  "Cl*”  ion,  for  which  only  one  broad  shell  of 
solvation  is  observed.  There  are  two  peaks  in  Cl'-hydrogen  rdf,  corresponding  to  two  different 
positions  of  hydrogens:  the  first  peak  is  due  to  the  hydrogens  that  are  oriented  towards  the 
ion,  the  second  is  due  to  hydrogens  that  are  pointing  away  from  the  ion.  In  the  case  of  *C1*”- 
hydrogen  rdf  the  picture  is  quite  different*  only  one  broad  peak  is  present  which  has  a 
shoulder  starting  at  ~5.7  A.  It  is  very  interesting  that  the  plot  of  coordination  numbers 
indicates  that  there  are  - 10  hydrogen  atnnm  beyond  thin  distance.  These  hydrogens  may  be 
considered  to  be  external  to  the  cluster  and  may  be  available  for  bonding  in  the  type  of 
titration  experiments  performed  by  Castleman’s  group  [6].  Note  that  in  the  idealized 
structure  proposed  by  Castleman  et,  aL  [6,8]  there  are  also  10  hydrogens  available  for  bonding. 
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Also  notice  here  that  the  peaks  of  "Cl+"-0  rdf  and  "Cl*"-H  rdf  appear  at  approximately  the 
mm*  distance,  although  one  normally  expects  to  see  the  cation-oxygen  peak  to  appear  at 
shorter  distances  than  cation-hydrogen  peak. 

Additional  information  on  the  structure  of  the  clusters  and  the  solvation  pattern 
around  the  ion  ««*  be  obtained  from  the  angular  distribution  function  (adf),  which  is 
displayed  in  Fig.  4.  In  this  figure  we  show  the  distribution  for  the  angles  between  the  vector 
connecting  the  ion  and  the  water  oxygen  and  the  vector  connecting  the  ion  with  the  center 
of  rnmut  of  water  molecules  that  are  the  nearest  neighbors  of  the  ion.  Thus  for  Cl*  ion  only 
waters  that  have  their  oxygens  closer  than  3.9  A  from  the  ion  are  taken  into  account  For 
the  "Cl**  ion  we  considered  all  the  water  molecules  that  have  their  oxygens  below  5.7  A, 
where,  as  we  mentioned  above,  a  beginning  of  a  slight  shoulder  in  the  water-ion  rdf  can  be 
observed.  As  we  can  see  from  Fig.  4,  the  behavior  of  the  adf  is  different  for  the  Cl'  and  "Cl*” 
ions.  The  adf  for  Cl'  ion  clearly  shows  that  the  ion  is  solvated  on  one  side,  while  water 
molecules  can  be  found  at  any  angle  around  the  "Cl*"  ion.  Note  that  even  for  "Cl*"  ion  the 
solvation  is  not  homogeneous,  pointing  out  that  also  in  this  case  the  ion  may  be  found  near 
the  surface  of  the  cluster. 

The  analysis  of  the  total  energy  of  the  clusters  and  its  components  can  help  us  to 
understand  the  structural  and  dynamical  behavior  of  the  dusters.  The  total  potential  energy 
of  the  duster  can  be  written  as  a  sum  of  two  components:  water-ion  and  water-water 
interaction  energies,  ie., 

Urn  -  <9> 

In  its  turn,  the  water-water  interaction  energy  can  be  further  divided  in  to  the  electrostatic 
(lPMrtl),  Lennard-J ones  (Uw‘wu)  and  polarization  (IPM’poJ)  components.  In  addition  to  these 
components  the  water-ion  interaction  energy  also  indudes  3-body  term  (IP’'4  M).  hi  Table  Ila 
we  display  the  time  averaged  values  of  all  the  components  for  the  water-ion  and  water-water 
interaction  energies.  In  Table  lib  we  present  the  total  values  of  the  electrostatic,  Lennard- 
Jones,  polarization  and  3-body  energies  and  finally  in  Table  lie  we  present  the  total  potential 
energies  and  the  two  components  of  this  energy  which  are  given  by  the  equation  (9). 

As  we  can  see  from  Table  II,  the  duster  with  Cl'  ion  is  the  most  energetically  stable 
one.  It  is  interesting  to  notice  that  the  contributions  of  water-water  and  water-ion  energies 
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into  the  total  energy  are  nearly  the  same  in  this  case.  Also  notice  the  rather  large 
contribution  of  the  polarization  energy.  Even  in  the  case  of  dewater  cluster  the 
polarization  energy  is  substantial,  mainly  due  to  the  polarizability  of  the  water  molecules. 
The  reduction  of  the  polarizability  on  the  C1‘  ion  results  in  a  slightly  weaker  ion-water 
interaction.  This  is  consistent  with  the  slight  change  in  the  shape  of  the  radial  distribution 
functions.  For  the  cluster  with  "Cl*"  ion  the  major  contribution  to  the  energy  comes  from  the 
water-water  interaction.  The  ion-water  interaction  in  this  case  is  weaker  compared  to  the  one 
observed  in  the  Cl*/ water  cluster,  in  spite  that  the  cluster  is  more  solvated.  The  strong  ion- 
water  interaction  observed  in  Cl*  water  duster  is  due  to  the  dose  proximity  of  water 
hydrogens  to  the  ion.  This  is  not  the  case  for  "Cl*"/water  duster,  where,  as  the  radial 
distributions  show,  oxygens  and  hydrogens  are  on  the  average  at  the  distance  away  from 
the  ion. 

To  get  a  better  understanding  of  the  relationship  between  the  duster  structure  and 
its  dynamics  we  hove  analysed  some  snapshots  of  duster  configurations.  Figures  la-lc  were 
used  as  guidelines  in  our  analysis.  As  Fig.  la  indicates,  the  Cl*  ion  is  expected  to  be  located 
on  the  surface  of  the  water  duster,  during  the  whole  trajectory.  This  was  confirmed  by  direct 
observations  of  duster  structures.  Some  representative  structures  of  Cl'fHjOja,  duster  are 
displayed  in  Fig.5.  Fig.  lb  indicates  that  dynamical  transitions  between  fully  solvated  and 
partially  solvated  structures  occur  in  the  case  of  "Cl*"  duster.  This  is  indeed  revealed  from 
the  snapshots,  which  are  displayed  in  Fig.  6.  While  the  ion  is  on  the  surface  of  the  duster  at 
t»332.9  ps  it  is  solvated  at  t- 377.1  ps  and  again  moves  to  the  surface  at  t* 422.4  ps.  The 
structure  at  t* 377.1  ps  is  dose  to  the  proposed  idealized  structure  which,  according  to 
Castleman  and  coworkers,  surrounds  Cs*  ion  [8].  For  comparison  we  present  the  two 
structures  in  Fig.7.  As  we  can  see  from  this  figure,  the  dynamical  structure  for  ”Cl*”/water 
duster  represents  a  torus  made  out  of  20  water  molecules  with  the  center  occupied  by  the 
ion.  The  idealized  structure  proposed  by  Castleman  et.aL  [6,8]  has  a  higher  symmetry  of  a 
pentagonal  dodecahedron. 

Minimum  Energy  Structures 

To  compare  the  energetics  of  the  idealized  "Cl*”(HsO)20  duster  structure  with  the 
energetics  of  the  other  stable  structures,  we  have  to  find  these  other  stable  structures.  For 
this  purpose  250  configurations,  selected  at  1  ps  interval  from  the  trajectories  of  the  dusters 
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with  Cl’  and  "Cl*"  ions  were  quenched  Figures  8  and  9.  represent  the  energy  components  of 
these  quenched  structures  which  correspond  to  different  local  minima.  A  salient  feature 
observed  from  these  figures  is  the  presence  of  a  broad  distribution  of  energies,  pointing  to  the 
existence  of  a  large  manifold  of  local  minimum  structures.  We  also  observe  from  the  figures 
that  the  distribution  of  water-ion  interaction  energy  for  the  quenched  Cl'/water  cluster  is 
broader  than  the  one  for  "CT  "/water  cluster,  indicating  that  the  clusters  with  the  Cl'  ion  are 
less  symmetrical.  This  is  consistent  with  the  conclusions  we  have  reached  previously. 

For  the  lowest  energy  configurations,  we  display  the  energies  and  their  components 
in  Tables  Ila-c  in  order  to  compare  them  with  the  average  energies  of  the  dynamical  dusters. 
As  we  can  see  from  the  Tables  and  Figure  8  the  water-water  interaction  energy  of  the  most 
energetically  stable  structure  of  Cl'/water  duster  falls  short  from  the  maximum  possible 
water-water  interaction  energy,  observed  in  quenched  structures,  but  it  is  not  very  far  from 
the  latter.  The  ion-water  interaction  energy  for  this  duster  is  right  in  the  middle  of  the 
observed  values  for  this  component  of  the  total  energy.  The  structure  corresponding  to  this 
minimum  energy  configuration  of  the  quenched  Cl'(H10)aD  duster  is  given  in  Fig.  10  a. 

For  the  most  stable  configuration  of  the  duster  with  "CT"  ion  the  water-water 
interaction  energy  is  at  its  lowest  value,  while  the  water-ion  energy  is  dose  to  the  highest 
value  observed  for  this  interaction.  Hie  structure  of  this  most  stable  configuration  of 
*Cl+"(H,O)J0  duster  is  displayed  in  Fig.  10b.  While  the  dynamical  configuration  has  a  shape 
of  a  torus,  we  can  see  that  the  most  energetically  stable  configuration  has  a  form  of  a  ’baseball 
glove11  with  the  ion  representing  the  "ball’. 

Hie  observed  minimum  energy  structures  from  Figures  10a  and  10b  and  a  comparison 
of  the  energetics  of  these  structures  with  the  data  from  Tables  II  indicate  that  the  stability 
of  the  quenched  structures  is  achieved  by  lowering  the  water-water  interaction  energy.  We 
also  observe  that  the  energy  of  the  idealized  structure  is  about  100  kJ/mol  higher  than  the 
energy  of  the  most  stable  quenched  structure. 

IV  CONCLUSIONS 

At  the  outset  of  the  work  we  raised  the  question  why  Cl'  ion  is  partially  solvated  by 
small  dusters  of  water  molecules.  To  check  if  that  was  due  to  a  large  ionic  polarizability  we 
have  investigated  the  structure  and  dynamics  of  dusters  where  the  polarizability  of  the  ion 
was  reduced.  To  see  the  effect  of  the  ionic  charge  we  investigated  the  structure  and  dynamics 


of  dusters  where  the  sign  of  the  ionic  charge  was  inverted.  From  our  simulations  we 
condude  that  the  reduction  in  the  ionic  polarizability  does  not  significantly  change  the 
structure  and  dynamics  of  Cl'CHjOljo  duster,  but  the  sign  inversion  of  the  ionic  charge 
produces  a  large  effect  While  the  Cl*  ion  prefers  to  be  situated  on  the  duster  surface,  the 
"Cl*"  ion  undergoes  dynamical  transitions  in  the  duster,  moving  from  its  surface  to  a  location 
whore  the  solvation  is  more  complete.  The  energetic  considerations  helped  us  to  understand 
why  Cl*  is  located  on  the  surface  of  the  duster.  By  being  on  the  surface  the  anion  permits 
the  creation  of  the  hydrogen  bonded  network  between  water  molecules  that  lowers  the  total 
energy  of  the  duster. 

When  the  solvation  of  the  ion  is  at  its  maximum  in  the  "Cr"(HjO)M  duster,  the 
duster  has  a  shape  of  a  torus.  This  structure  is  topologically  similar  to  the  pentagonal 
dodecahedron,  which  is  the  proposed  idealized  structure  for  the  Cs+(HtO)I0  duster. 
Nevertheless,  by  quenching  some  selected  configurations  obtained  from  the  dynamical 
simulations,  we  observed  that  the  minimum  energy  structure  has  a  shape  of  a  "baseball  glove", 
and  is  about  100  kJ/mol  lower  in  energy  than  the  idealized  structure.  In  the  present  work, 
our  nugor  goal  was  to  study  the  effect  of  the  ionic  charge  inversion  on  the  structure  and 
dynamics  of  the  ci'/water  dusters  and  not  the  structure  and  dynamics  of  Cs*/water  dusters. 
But,  since  "Cl*"  and  Cs*  ions  have  similar  sizes  and  polarizabilities,  we  believe  that  our 
conduaiona  reached  about  "Cl^'dijO)*,  duster  can  be  applied  as  well  to  a  Cs^dijO)^  duster. 

ft  is  interesting  to  compare  our  results  with  the  results  obtained  from  other  computer 
simulations  performed  on  inhomogeneous  aqueous  systems,  for  example  with  the  results  from 
the  simulations  that  study  the  structure  and  dynamics  of  the  ion  and  water  when  the  ion  is 
dose  to  the  water/vapor  interface  [18*201.  The  results  from  the  simulations  [18-20]  show  that 
water  molecules  at  the  interface  interact  with  the  anion  more  favorably  with  the  cation. 
We  have  reached  a  similar  condusion  in  our  simulations  performed  on  small  aqueous  dusters. 
But  contrary  to  the  observations  we  made  in  small  dusters,  it  was  observed  in  references  [18* 
20]  that  ions  keep  their  solvation  shell  at  the  water/vapor  interface,  ft  is  also  known  [18-20] 
that  the  ions  are  negatively  adsorbed  on  the  water/vapor  surface,  Le.  it  costs  energy  to  move 
the  ions  from  the  bulk  to  the  surface.  That  mown*  that  for  the  large  size  dusters  we  should 
observe  the  complete  solvation  of  the  anion.  The  question  that  remains  to  be  answered  is 
how  large  is  large? 
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FIGURE  CAPTIONS. 

FIGURE  L  The  time  dependence  of  the  distance  from  the  center  of  mass  to  the  ion  ((a)  Cl', 
(brer*  and  (c)  The  time  dependent  distance  between  the  center  of 
mass  of  the  water  molecules  and  the  whole  cluster  is  also  shown  and  it 
corresponds  to  the  curves  which  are  dose  to  zero. 

FIGURE  2.  Calculated  ion-oxygen  radial  distribution  functions  (bottom)  and  the  running 
coordination  numbers  (top).  The  solid  line  is  for  Cl'-O  rdf,  the  dotted  line  is  for 
"CT"-0  rdf  and  the  dashed  line  is  for  CT„-0  rdf. 

FIGURE  3.  Calculated  ion-hydrogen  radial  distribution  functions  (bottom)  and  the  running 
coordination  numbers  (top).  Solid  line  is  for  Cl*-H  rdf,  dotted  line  is  for  "Cr"-H 
rdf  and  the  dashed  line  is  for  Cl'^-H  rdf 

FIGURE  4.  Calculated  angular  distribution  functions  for  Cl'CHjOa,  (dotted  line)  and 
“Cl*"(HJO)I0  (solid  line)  dusters.  The  definition  of  the  angle  e  is  given  in  the 
text. 

FIGURE  5.  Snapshots  taken  from  a  trajectory  of  CI'CHjOIjq  duster.  Time  (in  ps)  at  which 
the  configuration  appeared  is  given  below  each  snapshot 

FIGURE  6.  Snapshots  taken  from  a  trajectory  of  "Q+"(HsO)20  duster. 

FIGURE  7.  a)  A  configuration  of  "C1+"(H,0)2o  duster  (at  t-377.1  ps).  b)  The  idealized 
structure  proposed  by  Castleman  etal[6,8]. 

FIGURE  8.  Ufc,  vs  Uw.w  (top)  and  U*,  vs  U,.w  (bottom)  for  the  quenched  Cl'CHjO)*, 
configurations. 

FIGURE  9.  U*.  vs  Uw.w  (top)  and  Uw  vs  U,.w  (bottom)  for  the  quenched  "Cr'GIjO)*, 

configurations. 
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FIGURE  10.  The  structures  corresponding  to  lowest  energy  configurations  obtained  in 
quenching  of  a)  Cl'dljO^  and  b)  "Cr'lHjO)^  clusters.  The  right  hand  side 
structure  in  (b)  is  the  side  view  (obtained  after  90°  rotation)  of  the  left  hand 
side  structure. 
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TABLE  CAPTIONS. 

TABLE  1.  Potential  energy  parameters,  charges  and  polarizabilities  used  in  the 
SPCE/POL  model. 

TABLE  2  a)  Components  of  the  potential  energy  (in  kJ /mol)  for  different  X(H20)2o  clusters. 

(The  subscripts  "rp",  •min"  and  "id*  are  used  to  represent  the  reduced 
polarizability  of  the  C1‘  ion,  the  lowest  minimum  energy  structure  found  during 
the  quenching  of  250  samples  and  the  idealized  structure  according  to 
Castleman  at  aL[6,8]  respectively.) 

b)  Energy  components  as  in  Equation  (1). 

c)  Total  potential  energies  and  their  major  components  as  in  Equation  (9). 
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TABLE  I. 


A*10*5  (Kcal  A12/mol)  C(kcal  A6/mol) 


0-0  6.37  627.95 


0-Ion  37.90  1373.84 


q(e)  <*(A3) 


-0.73  0.465 


H  0.365  0.135 


Cl*  -1.0  3.25 


"Cl*"  1.0  3.25 


Cl*,.,,  -1.0  0.0025 


B*104(kcal/mol) 

BtA*1) 

Y^*1) 


-  0.80 
m  2.25 
-  0.25 


16 


TABLE  II. 


(a) 


Cl* 

■  cr* 

Cl* 

fp 

Cl'rfn 

"cr%ln 

"Cl*",. 

U--.1 

-593.4 

-644.1 

-562.1 

-776.4 

-804.5 

-772.0 

205.8 

227.9 

188.7 

318.5 

284.8 

223.6 

-74.5 

-125.9 

-83.1 

-120.8 

-173.4 

-120.9 

0'-“., 

-256.6 

-372.8 

-376.0 

-250.8 

-154.0 

°'*u 

53.1 

11.8 

37.4 

64.0 

8.1 

-5.4 

°'-V 

-132.1 

-52.0 

-78.0 

-162.9 

-55.0 

-56.2 

rjf-M 

u  3-b 

8.0 

5.0 

7.2 

8.9 

6.9 

5.3 

(b) 


Cl* 

■cr» 

Cl* 

C*  *ln 

"Cl*".,„ 

"Cl*" 

id 

-954.1 

-900.7 

-934.9 

-1152.4 

-1055.3 

-926.0 

UU 

258.9 

239.7 

226.1 

382.5 

292.9 

218.2 

V 

-206.6 

-177 . 9 

-161.1 

-283.7 

-228.4 

-177.1 

°3-b 

09 

• 

o 

B 

7.2 

8.9 

6.9 

5.3 

(c) 


Cl* 

"Cl*" 

Cl* 

Cl«. 

"Cl*"id 

-462.1 

-542.1 

-456.5 

-578.7 

-693.1 

-669.3 

-431.7 

-291.8 

-406.2 

-466.0 

-290.8 

-210.3 

"tot 

-893.8 

-833.8 

-862.9 

-1044.7 

-983.9 

-879.6 

inYrnainom^FicD^om^tnai^  oc\j 


(V)  aDue^sia 
© 


Time  (ps) 


!•!»*»* 


40- 


350.0  380.0  430.0 


096- 


xeioin 


le^o^n 


1010 


-006- 


1000 
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